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Low temperature dynamics of the S = ^ Heisenberg chain is studied via a simple ansatz gener- 
alizing the conformal mapping and analytic continuation procedures to correlation functions with 
multiplicative logarithmic factors. Closed form expressions for the dynamic susceptibility and the 
NMR relaxation rates 1/Ti and I/T2G are obtained, and are argued to improve the agreement with 
recent experiments. Scaling in q/T and uj/T are violated due to these logarithmic terms. Numerical 
results show that the logarithmic corrections are very robust. While not yet in the asymptotic low 
temperature regime, they provide striking qualitative confirmation of the theoretical results. 

PACS: 75.10.Jm, 75.40.Gb, 75.50.Ee, 76.60.-k 



O 



> 

in 

o 
o 

a 

O 

o 



X 
S3 



In recent years there has been much interest in quan- 
tum critical phenomena in spin-models and real mate- 
rials. In path integral formulations, the inverse tem- 
perature /3 acts as a finite size for the imaginary-time 
fluctuations, driving the system away from the T = 
quantum critical point. The resulting behavior can be 
described in a manner analogous to finite-size scaling. 
In two-dimensions, quantum critical points are rare, but 
their relevance to real materials is enhanced by the fact 
that they also control the finite temperature properties 
of weakly-ordered or weakly-gapped systems . 

In contrast to 2D, ID quantum antiferromagncts with 
continuous symmetry are generically critical at T — 0. 
Thus one expects to find many examples of quantum 
critical phenomena in quasi-one dimensional materials. 
Recently, quantum critical behavior of spin-chains have 
been studied by neutron scattering and NMR experi- 
ments Q. 

From a theoretical point of view, the development of 
conformal field theory provides a powerful machinery to 
study the finite temperature correlation functions at a 
T = critical point. Assuming a simple power-law be- 
havior for the T = correlation functions, various au- 
thors have obtained scaling forms for the dynamic struc- 
ture factor at low temperatures However, it is 
well known that the T = spin-spin correlations of the 
S = 1/2 chain, have multiplicative logarithmic factors 
due to the presence of marginally irrelevant operators 
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Here c is the spin wave velocity, and D and tq are nonuni- 
versal constants. These logarithmic factors also appear 
in the two-spinon contribution to the dynamic structure 
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In this paper, we explore the effects of logarithmic fac- 
tors on the finite temperature dynamics of the spin-half 
Heisenberg chain. We begin by considering the T — 
equal-time correlations of a finite system. It was pro- 
posed in Ref. that, in presence of logarithms, the cor- 
relations for a system of size L should have a generalized 
finite-size scaling form: 
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Here the universal function X{x/L) is given by X{x) = 
i sin(7ra;) [|l^. This proposition was shown to work well 
for the isotropic Heisenberg chain with cr = 1/2 

In the absence of logarithmic factors, correlations at 
finite T are obtained from a conformal mapping of the 
complex plane z = x + icr at T = to a strip infinite in 
the x-direction and of width c/T in the r-direction (r is 
imaginary time). In close analogy with T = finite-size 
scaling, this amounts to the substitution 
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We now assume that this mapping can also be performed 
in presence of marginal interactions. Then finite T spin 
correlations are given by the ansatz (H) with scaling func- 
tion (|^). This can be further approximated as: 
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where Tq = v^^rc/ro, and an effective temperature 
dependent scaling dimension appears [pd|,p[ 
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Equation (^) is valid for a; <C f In ^ (see expression for 
^ below), and thus allows one to study correlations be- 
low and above the correlation length. An immediate 
consequence of this ansatz is that the correlation length 
aquires a logarithmic temperature dependence, in agree- 
ment with thermal Bethe ansatz calculations [|l2| : 
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We now explore further consequences of the scaling 
ansatz. First, the static structure factor is found to be 



S{q) = 2^^+'^^D In 



To 
T 



1/2 



r(i-4A) 



xRe 



r(2A- 



'2ttT) 



r(i-2A 



''2itT I 



(7) 



where q is measured from the antiferromagnetic vector 
TT. Note that the entire q-dependence of S{q) is due to 
the 1/lnT corrections to the T = value of A = 1/4. 
Eq. (0) implies that S{ci)/S{Q) is no longer a universal 
function of cq/T . 

Performing Fourier transformation and analytic con- 
tinuation to real frequencies |§,||,D, we obtain the stag- 
gered susceptibility 
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This expression also lacks universality due to the T- 
dependence of A. 

Next, we discuss numerical results for the spin-half 
chain obtained using a "stochastic series expansion" 
Quantum Monte Carlo (QMC) method (for systems 
with up to 1024 spins), and conventional high tempera- 
ture expansions (HTE) . Most results from the two meth- 
ods agree down to T/J =1/8. Below that temperature, 
we rely on QMC data alone. 

We begin with the w = susceptibility, shown in Fig. ||. 
The ratio x((7, 0)/x(0, 0) appears to converge towards a 
scaling form as the temperature is lowered, but even at 
/3 = 32 it is far from the universal scaling function ex- 
pected in the absence of logarithms [Q. In the range 
l/4>r>l/8 the numerical results have high accuracy, 
and QMC and HTE data agree very well. The deviations 
from scaling are clearly systematic, and well described by 
Eq. (^) with Tq = 4.5. Note that the parameter Tq should 
be considered an effective one. As the study of logarith- 
mic corrections to the uniform susceptibility shows [ p^ , 
the true value of Tq maybe reached only at T < 0.01. 



Data for S{q) show substantial g-dependence, in dis- 
agreement with A = 1/4 scaling predictions. However, 
the results are not well explained by Eq. (Q) either. A 
possible reason is that S{q) is dominated by contributions 
(divergent at T = 0) from short distances, where our 
asymptotic expression (^) breaks down. It is, thus, better 
to compare the equal-time real-space spin correlations, 
S{x), with the theoretical expressions. It is well known 
that the correlation function in addition to the domi- 
nant staggered piece has a uniform contribution, given 
~( 2csinhf7rTa;/c) )^ ^"^^^^ T' Q • ^ IS appropriate to 
subtract this from the numerical data before comparing 
with the scaling theory. As shown in Fig. |^, our results 
for iS'(a;) agree very well with Eq. (^, with Tq = 4.5 and 
D = 0.075. The inset shows a comparison of the ra- 
tio of correlation functions at two temperatures. With 
To fixed from the susceptibility data, this parameter-free 
agreement is quite striking. Deviation of the theoreti- 
cal results at short distances is also apparent and is the 
reason that S{q) cannot be explained. The theoretical re- 
sults also imply S'(O) ~ (In/?)^/^ and x(0,0) - /3(ln/3)i/2 
as T 0, in agreement with numerical data [l^ . 

From Eq. (||) we can calculate the NMR relaxation 
rates 
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Here the integrals Ti(A) = dxj^r^^^^ auu ^ 

^ /o°° I r(i^A-I) 1'^ have weak temperature dependences. 
In deriving Eq. (0), we have kept only the scaling part 
and dropped the term coming from self interactions as 
it is down by a factor T(ln^)^. The latter is, in any 
case, not correctly accounted for by the scaling theory. 
We find that Eq. (||) shows weaker than ^JIyiTq/T vari- 
ations with T. This result is in qualitative agreement 
with recent measurement of 1/Ti in Sr2Cu03 by Taki- 
gawa et al. js). Fig. || shows the ratio T^gI \fTT\. In the 
T ^ limit our expressions (||) and ([l0| ) coincide with 
those of Sachdev |lj]. However, we find that T = limit 
of T^gI^/TT^i is approached with infinite slope, similar 
to the behavior of the uniform susceptibility jl^. The 
behavior is consistent with the slow rise of this quantity 
seen for SraCuOg around T = J/10 |]. 

Numerical results for x" are obtained from QMC data 
using the maximum-entropy (max-ent) method p], and 
from HTE via the recursion method In Fig.^, data 
are presented for the ratio with the full T^g and with only 
the scaling part, where the self-term is not subtracted 
p8[ . The two should converge in the scaling limit and 
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the latter should be compared with the theoretical re- 
sult. Note that QMC and HTE results agree completely 
for T2G', the deviations arise entirely from the analytic 
continuation needed to get l/Ti, which is more uncer- 
tain for QMC data (details of this point will be discussed 
elsewhere Q). The difference between the curves based 
on the full T2G and the scaling part only of T2G shows 
that the results are not yet in the scaling limit. However, 
the theoretical results are supported by the convergence 
of the more accurate (in the temperature regime shown) 
HTE data to the predicted form. The presence of non- 
asymptotic contributions in the full T2G and the apparent 
tendency of QMC -I- max-ent to over-estimate 1/Ti ex- 
plain the discrepancy in the numerical result previously 
reported for T2G/VTT1 jl9|. 

We also note that the experimental 1/Ti was found 
to be about 30% lower than the numerical result at 
T = 300K [|jl9|, which is now also reconciled. Together 
with the good agreement found previously for I/T2, with- 
out adjustable parameters ||3|,|l9|, the spin- half chain in- 
deed very well describes the low-frequency dynamics of 
SraCuOs. 

The frequency-dependent quantities also do not show 
universality in the scaled variable ui/T. For example, 
the imaginary part of the q — and local (xi^) — 
1^00 ^xilj^)) susceptibility are given, respectively, by 



to explore the scaling forms at larger oj/T. There are 
preliminary reports of measurements of these quantities 
by neutron scattering |^ . It would be useful to compare 
them with our results. 

To conclude, we have studied the effects of logarith- 
mic corrections on the finite temperature dynamic spin- 
correlations of the spin- half chain. Analytical expressions 
are developed for x(9, w) by a generalized finite-size scal- 
ing ansatz. The ansatz ties together previous results, in- 
cluding logarithmic corrections to the correlation length, 
and implies a temperature dependent effective scaling di- 
mension. Expressions obtained for the NMR relaxation 
rates are argued to improve the agreement with experi- 
mental data for Sr2Cu03 Numerical results, although 
not in the asymptotic low temperature regime, confirm 
various theoretical expressions including violation of scal- 
ing in the variables cq/T and oj/T. We expect these ef- 
fects to diminish and scaling to be restored if the second- 
neighbor interactions are tuned to the point where the 
marginal interaction is absent We hope our calcu- 

lations will provide further motivation for neutron scat- 
tering studies of quasi-one dimensional spin systems. 
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The a;-dependence in Eq. ( p^ is due to the T-dependence 
of the scaling dimension. The a;-indcpendence and di- 
vergence at A = 1/4 is an artifact of the approxima- 
tions employed, and should be removed by more accu- 
rate treatment of the short-distance cut-off. Neverthe- 
less, fixing the scaling dimension at 1/4 implies a much 
weaker T— and tu— dependence than predicted by (12). 
Note also that we predict x"(0,^) ^ (In^y/^ and 
x"(a;) ~ (In^)"^/^ at low temperatures. 

QMC -I- max-ent data for the q = susceptibility, as 
well as the results of Eq. (|ll|) , are shown in Fig. ^, where 
the value of Tq is from the fit of the static susceptibility. 
At low frequencies the theory and the data agree and also 
appear to scale. At higher frequencies, there is no scaling 
and the deviations are qualitatively similar in that the 
lower temperature data is higher at higher values of uj/T. 
The numerical data are not at low enough temperatures 
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FIG. 1. The static susceptibility normalized to its 
g = value. Symbols represent numerical data from 
high-temperature expansions (HTE) and Monte Carlo sim- 
ulations (MC). Solid lines are predictions of Eq. (^ with 
To = 4.5, and the dashed line shows the universal scaling 
function with A = 1/4 [|l^. The inset shows the QMC data 
over a larger range of f3 and cq/T together with the universal 
scaling function. 



FIG. 2. Comparison of QMC data for real-space correlation 
functions ( symbols ) and Eq. (^) ( solid lines ) with To — 4.5 
and D = 0.075. The inset shows the ratio of the equal-time 
correlation functions at /3 = 16 and 32 compared with Eq. (^) 
and the expression with A = 1/4 [Q. 



FIG. 3. The ratio T2g/VtTi versus T from Eqs. (|) and 
(0). The T = limit of the ratio is 1.68. The inset shows 
a linear variation with 1/ln {[3). QMC and HTE data for the 
ratio with full T2G and with only the scaling part included 
are shown by the symbols. 



FIG. 4. Imaginary part of the antiferromagnetic suscepti- 
bility. Lines represent Eq. ( pi] ) with parameters as in Fig. |^ 
and the symbols represent the QMC -I- max-ent data. 
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